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We find a matter-wave solution of Bose-Einstein condensates trapped in a harmonic-oscillator 
potential and subjected to a homogeneous gravitational field, by means of the extended tanh- 
function method. This solution has as special cases the bright and dark solitons. We investigate 
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without changing the strengths of the magnetic and gravitational fields. 
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INTRODUCTION 



Bose-Einstein condensation is a very general physical 
phenomenon which takes place in the systems of bosonic 
atoms at ultralow temperatures, as well as in optical wave 
systems It appears in the fields of condensed mat- 
ter, atomic and elementary particle physics and astro- 
physics Q. The majority of important features of the 
condensation in such diverse systems can be captured 
by the Gross-Pitaevskii (GP) equation, which is a vari- 
ant of the nonlinear Schrodinger equation with a trap 
potential. Due to the nonlinearities arising from the in- 
teratomic interactions and due to the presence of a con- 
fining potential in the GP equation, many studies have 
been performed either by solving the corresponding GP 
equation numerically or by using perturbative methods 
(see [3|, |4{ and references therein) . However the construc- 
tion of exact solutions of the GP equation can help to ad- 
vance our understanding of the various physical phenom- 
ena governed by this nonlinear equation. For instance, 
the exact analytical solutions can contribute to select 
the experimental parameters, to analyze the stability of 
Bose-Einstein condensates (BECs) and to check numeri- 
cal analysis of this nonlinear equation [5[- Therefore, the 
construction of the exact solution of the GP equation is 
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one of the most relevant challenges to take up. To this 
end, many works have proposed several methods to ex- 
actly solve the GP equation, such as the Darboux trans- 
formation method, the hyperbolic-function method, the 
elliptic-functions method, the inverse-scattering method, 
the generalized (G /G)- expansion method, the Hirota bi- 
linear method, the Painleve analysis, the Lie symme- 
try, the tanh-function method and the extended tanh- 
function method (see for instance |6j and references 
therein) . Finding the nontrivial seed solutions that mim- 
ics the sought solutions of the GP equation always re- 
mains a hard task. 

Recently, there have been some experimental reports 
on solitons in the quasi one-dimensional (ID) system 
which have been observed by magnetically tuning the 
interparticle interaction from repulsive to attractive. Ex- 
perimentally, for the realization of ID systems, one 
should make the radial frequency much larger than the 
axial frequency and strongly confine the radial motion. 
The state in this cylindrical harmonic-oscillator trap rep- 
resents a nonlinear matter wave. 

The formation and propagation of matter waves, such 
as dark solitons, bright solitons, four- wave mixing and 
moving or stationary gap solitons [3] are among the most 
interesting dynamical features in Bose-Einstein conden- 
sation. Of course, such dynamics depends on the types 
of interactions in which the system is subjected to. For 
atoms in the nK-mK temperature regime, the effect of 
the Earth's gravitational field is by no means negligible 
especially in the case of magnetic trapping. The gravita- 
tional field plays a subtle role in the topological formation 
of stable vortices within the field reverse time in BEC ex- 
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periments with heavier atoms like 87 Rb It has been 
shown that the gravitational field can change the propa- 
gation trail of the bright soliton trains without changing 
the peak and width of the soliton in the parabolic back- 
ground [5|, [9j. Moreover, the presence of a homogeneous 
gravitational field can decrease the condensation temper- 
ature of Bose gases [lOj . 

BECs also appear in the astrophysical context. For 
example, one has recently argued that dark matter 
an unknown component which corresponds to approxi- 
mately 25 % of the energy content of the Universe, is 
currently under the form of a BEC [12]. In this approach, 
bosonic dark matter particles underwent a phase transi- 
tion forming a BEC at some point during the evolution of 
the Universe. As a consequence, one opens the possibil- 
ity of admitting the existence of BECs subjected to many 
gravitational effects. This seems to be a very promising 
line of research for the next years. 

In this paper, we study the dynamics of BECs 
in a magnetic trap in the presence of gravity. We 
construct an exact solution of the GP equation which 
has as special cases the well-known bright and dark 
solitons. The work is organized as follows. In Sec. 
HIl we present the model. Then in Sec. IIII| by uti- 
lizing the extended tanh-function method, we find 
the various soliton solutions. In Sec. IIVI we examine 
the dynamical properties of these solutions. Finally in 
Sec. El we summarize our results and conclude the work. 



II. THE MODEL 

As is well known, in the mean field approximation, the 
full dynamics of a BEC in a trap potential V(r) satisfies 
the time-dependent Gross-Pitaevskii (GP) equation fl3l — 

m 

ifitft = -t^V 2 ^ + V(r)tf + £o|^| 2 ^ (1) 
2m 

where go = 47r ^° a , with m being the atomic mass and a s 
the 5-wave scattering length which can be either pos- 
itive (the case of 87 Rb atoms with a s = 5.45 ± 0.26 
nm and then repulsive interactions) [16| or negative (the 
case of 7 Li atoms with a s — —1.45 nm and then attrac- 
tive interactions) [l5j. It has been shown that, for al- 
kali atoms at least, the effect of gravity is non- negligible 
We recall that in the usual experimental traps, 
the atomic clouds are confined with the help of laser or 
magnetic trapping. For alkali atoms such as rubidium, 
which is the most used boson for experiments on cold 
atoms, most of the existing traps can be well approxi- 
mated by a three-dimensional harmonic oscillator [14j . 
We consider a trapped Bose gas of non-relativistic par- 
ticles immersed in a Newtonian gravitational field. In 
this case, the external confining potential must take into 
account the contribution of both the magnetic trapping 
field and the gravitational field. It is given by the sum of 



the quadratic and gravitational potentials generated by 
these respective fields [l(| l2Q| : 

V(r) = j (uj 2 x x 2 + uj 2 y y 2 + uj 2 z 2 ) + mgz, (2) 

where lji, with I = x,y,z, denote the frequencies of the 
harmonic oscillator along the coordinate axes. The pa- 
rameter g, which represents the acceleration of gravity, 
is taken as constant since the gravitational field is homo- 
geneous. The gravitational potential represented by the 
last term of Eq^ (|2]) is also present in the case of the 
tilted trap HEElI- 

In general, the harmonic-oscillator potential repre- 
sented by the first term in Eq. ((2]) is anisotropic, i.e., 
the trapping frequencies uji are all different. For a trap 
that is axially symmetric along the z- axis, we have 
w x = = u±. In such a case, uo z is referred to as the lon- 
gitudinal frequency (along the axial direction) while uo±_ 
is the radial frequency of the anisotropic harmonic trap. 
Then Eq. (pQ) reduces to the following three-dimensional 
GP equation: 

h 2 1 

ifi*t = -— V 2 ^+[-m(cjiri+cj^ 2 )+m^]^+^ |^| 2 ^, 

m (3) 

where r± = x 2 + y 2 is the radial distance. The radial 
motion can be strongly confined by making the radial fre- 
quency much larger than the axial frequency, i.e. uj± ^> 
uo z . In this case the condensate is cigar-shaped, and ow- 
ing to that [22|, one can take \I/(r, t) = (/>o(r±)ip(z,t), 

where 0o = exp( — ^-) is the ground state of the 

radial problem, with a± = ^/h/muj^. Then multiplying 
both sides of the GP equation (j3j) by 0q and integrating 
over the transverse variable r_\_ we obtain a quasi-one- 
dimensional GP equation in the form: 

h 2 1 

ihi/; t = —^l>zz + (-mu^z 2 + m 9 z )^ + v\i/j\ 2 i/j. (4) 

Thus, the GP equation (j4]) describes the dynamics of 
trapped quasi- ID cigar-shaped Bose gases at the mean- 
field level. In this equation, the strength of the atom- 
atom interaction becomes v = 2huo±a s . The s- wave scat- 
tering length a s can be managed through the Feshbach 
resonance technique [23|. Additionally, the effect of the 
acceleration of gravity g in the BEC experiment can be 
tiny varied or even cancelled in drop tower experiments 
[241 ]. In the astrophysical context, we know that BEC 
might exist (in a speculative scenario) and can be sub- 
jected to very high gravitational fields (close to black 
holes, for example). This suggests the possibility to vary 
the gravitational field. Hence one has some freedom in 
choosing the physical parameters of the system. In what 
follows, we set = c, \muj 2 = ha, mg = /iA, and 
v = hv{t). 

We follow the ideas of Refs. 0, H| and seek the exact 
solitonic solutions of Eq. ((4]) within the extended tanh- 
function method. We therefore introduce the auxiliary 
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equation 



>2 



C + C\U + c 2 u 2 + c 3 ^ 3 + c 4 ^ 4 , 



(5) 

where ^' = C = C( z ^) and Co, ci, C2, C3, C4 are some 
real constants. Let a be an arbitrary real constant. One 
may verify that Eq. (j5j) admits for Co = 0, c\ = 0, C2 > 0, 
and C3,C4 = arbitrary, the following solution: 



«(C) 



4c 2 



-2c 3 + (c§ - 4c 2 c 4 )e' 5 v / ^^- a ) + e-^v^lC-a) ' 



(6) 



where S = ±1. It is worth noting that this solution 
presents some singularities for c 2 — 4C2C4 < 0. 

We let C(z : t) = p{t)z-\-q{t) and use the transformation 

^(*, t) = [/(*) + /i(t)^(C)]e i ^ (t) ^ +fc(t ^ +Q(t)] . (7) 

The parameter a is called the steepness of the harmonic 
trapping potential, in reference to the potential energy of 
a perfectly elastic spring. f3(t) is twice the chirp and Q(t) 
is the linear phase. Next, we substitute Eqs. (j5]) and ([7]) 
into Eq. ^ , and set the real and imaginary parts of the 
resulting equation to zero. Then collecting coefficients of 
powers of z n u n (n = 0, 1; n' = 0, 1, 2, 3, 4, 5) and setting 
each of them to zero, we obtain the following set of over- 
determined ordinary differential equations: 

/' + 2c(3f = 0, ti + 2c(3h = 0, 
h(p' + 4cf3p) = 0, h(q' + 2ckp) = 0, 
h{k' + 4c/3/c + A) = 0, f(k' + 4c/3fc + A) = 0, 
/(/?' -f 4c/3 2 + a) = 0, h(p + 4c /3 2 + a) = 0, 



(8) 



h(2cp 2 c 4 - vh 2 ) = 0, -h(cp 2 c 3 - 2vfh) = 0, 



2f(Q f + ck 2 + vf 2 ) 
h(VL' + ck 2 + Zvf 2 



- c hp 2 ci = 0, 
cp 2 c 2 ) = 0. 



In this set of equations, the prime stands for the deriva- 
tive with respect to the time t. Since c\ = is the con- 
straint on the solution (j6j), the eleventh equation of the 
set of equations J8J) may read f(Q' + ck 2 + ^/ 2 ) = 0. 
Applying different constraints between the parameters 
in Eq. (|6j), the above equations may yield to various 
solutions. Thus to solve Eq. (|5J), we consider two situ- 
ations that have interesting implications on the matter 
wave solution, namely when C3 = and C3 ^ 0. 



III. THE SOLUTIONS 

A. The bright solitons 

For C3 = 0, we get f(t) = 0, and then we may reduce 
the above set of equations J5J) to 



ti + 2cfih = 0, p' + 4c/3p = 0, g' 
fc' + 4c/3fc + A = 0, P' + 4c/3 2 



2ckp - 
a = 0, 



0, 



(9) 



ft' + c/c 2 



cp 2 c 2 



In the case of time-independent magnetic and gravita- 
tional fields, say a and A are constants, we may choose 
the parameter f3(t) in such a way that the above set is 
easily solvable. We consider the following three exam- 
ples. 

(1) When /3(t) is constant, solving the above set of 
equations yields: 



k(t) = 

<?(*) = 
fi(t) = 



~,p(t)=p £(t)-\h(t) = h £(t)- 1 / 2 , 
4c 



4c/3 



+ h£(t)-\ v(i) = u £(t)-\ £(t) = e 



-4c(3t 



A 2 



16c/? 2 



C2P0 ~ fc| Sc 0t _ Mi_ 



8/3 



8cf3 2 



lit)- 1 

(10) 



In this set of equations, ^0 and fto = ^^c^p 2 ,) 1 / 2 are 
the initial values of z/(t) and h(i), respectively. The pa- 
rameter po is a nonzero real constant while fci, gi and 
fti are some arbitrary real constants. These results im- 
plicitly assume vqc± > and particularly a < 0, in view 
of the above expressions of f3 and feo- Thus the solution 
in this case corresponds to expulsive magnetic trapping 
potential. Furthermore, it demands C4 ^ in order for 
the amplitude coefficient h to be nonzero. 

(2) When f3(t) = -yff c tan(2 v /cat), solving the set of 
equations (|9]), we come to the following results 

p(t)=p £(t)-\h(t) = ho£(ty 1/2 , 
A 



k(t) 



2\/ca 



tan(2vcat) + 



-ki\ — sm(2^/cat) + ^~ 
V a 2a 



P(i) + qi, 



~, x A 2 A 2 - 4ca(c 2 pn - &?) /rt / — x 

=-— t 5l_^o IL tan(2vcat) 

v ; 4a SaVca v v 7 



0, 2cp 2 c 4 - vh 2 



0. 



(11) 



=z/ ^(t) _1 , ^(t) = |cos(2vc^t)|. 

These results also assume ^0^4 > 0, with C4 ^ and 
especially a > 0. Hence the solution in this second case 
corresponds to attractive magnetic trapping potential. It 
signals singularities at any t = , with n being any 

positive integer. In such case we consider only safe times 
defined by t < to prevent singularities. This as- 

sumption does not alter the validity of our approach since 
the lifetime of a condensate is small in general. Moreover 
most of experiments are done with weak magnetic field, 
i.e., with small values of the parameter \a\. In this situ- 
ation, the cut-off time t\ = is long enough to allow 
the observation of the matter-wave propagation in the 
condensate. 
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(3) When /3(t) = ^/~^ tanh(2 v / 3 ca t), solving the 
set of equations ([9]), we obtain the following results 

p{t) =pd{t)~\ h(t) = hoeit)- 1 ' 2 , 
A 



k(t) 
«(*) 



tanh(2\/-cat) + M(t) -1 
p(*) +9i, 



\2 

nrt) =— t 

4a 



2\/— co; 

a 2a 
A 2 - 4ca(c 2 pl - fc 2 ) 



(12) 



Afci 
~2a 



£(t) 



SayJ—ca 



tanh(2V— cat) 



v(t) =i^(*)~\ ^) = cosh(2v^cat). 

Likewise these results also assume ^o c 4 > 0, with C4 7^ 
but especially a < 0. Hence the solution in this third case 
corresponds to expulsive magnetic trapping potential. 

Considering the above results, the general solution of 
equation (j4| for C3 = may be given by 



i/>(z,i) = h(t)u(£) e 1 



((3z 2 +kz+Q) 



(13) 



where u(Q is the solution of Eq. ([5]) given in Eq. (|6]). 
In order to obtain a soliton solution from Eqs. ((6]) and 
(Q3]), we need to avoid singularities in Eq. (|6j). Then, we 
assume in Eq. (J6j) the constraint C4 < 0. In this case, 
taking C3 = as demanded above, one can readily obtain 



u(() = u sech[^(( - Co)], 



(14) 



where uq 



and Co 



= ln(— 4C2C4). In this 



C4 ^ W ~ 2^c 2 

case Eq. (fT3j) describes a bright soliton. This solution is 
displayed in Fig. [TJ It should be noted that since C4 < 0, 
we have negative v$ and then, this bright soliton is meant 
for media with attractive interparticle interactions. 



B. The dark solitons 

Another case that may be of interest in obtaining soli- 
ton solutions is when cs 7^ 0. We may have f(t) 7^ 0, and 
reduce the set of equations (|5]) to 

ti + 2c/3h = 0, p f + 4c/3p = 0, q' + 2c/cp = 0, 

fc' + 4c/3fc + A = 0, /3' + 4c/3 2 + a = 0, 

2cp 2 c 4 - z^ 2 = 0, cp 2 c 3 - 2vfh = 0, f' + 2cf3f = 0, 



tt' + c/c 2 + -cp 2 c 2 = 0, cp 2 c 2 - 2vf 



0. 



From Eq. ([T5]) . it may be shown that, when C3 7^ 
the parameters in Eq. (j6j) must fulfill the constraint 
c 2 = 4C2C4. In the case where a and A are constants, we 
choose the parameter j3(t) as in the previous section to 
easily solve the above set of equations. 






FIG. 1: (Color online) Plot of the solution (p) for (a) 
= \F¥d ( b ) PV) = y^tanh(2 v / ^t), and (c) 
Pit) = —\p^. tan(2y / ca£). The parameters are a — ±0.005, 
po — 1.0, iso — 1.0, a = 0, c = 0.5, C2 = 1, C4 = z/o, &i — 1.0, 
<7i = 0, and A = 0.01. Then the cut-off time is approximately 
ti — 15.7. All quantities are dimensionless. 



(1) When /3(t) is constant, the solution of Eq. ([T5]) is: 



^,p(t)=p e-^ t , h(t) = h e- 2c ? t , 



f{t) = foe- 2 ^, uit) = u e- 
X 



-4c/3t 



k{t) 



4c(3 



-4c/3t 



'("=(-^ + l e "^) p(t)+ 

n(t) 



9i, 



16c/3< 



C2P0 + gfci -Sc/Bt 

+ 16/? 



Afci 
8^ 



" 4c ' 3t + fi 1 . 
(16) 



(15) 

In this set of equations, fo 



f^ho is the initial value 
of fit). The difference between Eq. ([T6]) and the cor- 
responding results in the previous section resides in the 
expressions for fit) and Q(t) which change, respectively, 
the amplitude and the phase of the matter wave. 

(2) When pit) = -J¥ tan(2 v /cat), solving the set of 
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equations (fT5]h we come to the following results 

P (t) = PQ e(t)-\ Kt) = m(*)- 1/2 , /(*) = M*r 1/2 , 

A 



fc(t) 



2 x [ca 



ta,n(2y/cat) + 



-k\ \ — sin(2 Vca t) + — 



A 2 A 2 + 2ca(c 2 p 2 ) + 2k?) , j— N 
=— t V — - tan(2-s/cat) 



z/(f) =zV(*) _1 , ^) = |cos(2\/cat)|. 



(17) 



We can prevent singularities by considering only the 
times preceding the cut-off time, i.e., we take t < t\. 
However this particular singularity can be avoided by 
changing the scattering length in a small time interval 
that contains each singular time t n . 

(3) When f3(t) = tanh(2v^ca t), solving the 

set of equations (fT5]h we obtain the following results 

p(t) = P0 £(t)-\ hit) = M(t) _1/2 , f(t) = f m~ 1/2 , 
A 



kit) 

nit) 



tanh(2 v /z cat) + 

pit) 



2^—ca 

C _0 . / — nr* + A 

a 2a 



A 2 A 2 + 2cajc 2 pj + 2k\) 

-—t . — tanh(2y— cat) 

^a Sa^—ca 



v (t) =i/ ^(t) _1 , £(t) = cosh(2\/^cat). 



(18) 



Considering the above results, the general solution of 
equation (j4]) for C3 7^ is given by 



i/>{z,t) = h(t) 



_C3_ 

4c 4 



■«(0 



i((3z 2 +kz+Q) 



(19) 



where ia(C) is the solution of Eq. (|5j) given in Eq. (|6j). 
The condition c 2 = 4C2C4 implicitly assumes C4 > since 
C2 is positive. In this case, we easily come to 



«(C) 



4c 2 



-2c 3 + exp[-(5^/5J (C - a)] ' 



(20) 



Eq. ([20]) is an interesting solution without singularity 
only when the nonzero coefficient C3 is negative. This 
solution, with Eq. ([T9]h corresponds to a darklike soliton. 
We portray this solution in Fig. [2j We mention that since 
C4 > 0, the parameter v$ is positive and consequently this 
dark soliton is meant for media that present repulsive 
interparticle interactions. 




FIG. 2: (Color online) Plot of the soluti on (H9l) for (a) 
P(t) = V^E> ( b ) PV) = V^^MV^S*), and (c) 
Pit) = — \p^. tan(2y / cat). The parameters are a = ±0.005, 
po — 1.0, — —1.0, a = 0, c = 0.5, c<i — 1, C4 = ^0, k\ — 1.0, 
gi = 0, A = 0.01, and 6=1. The cut-off time is roughly 
t\ = 15.7. All quantities are dimensionless. 



IV. THE DYNAMICAL PROPERTIES OF THE 
SOLUTIONS 

In order to investigate the dynamical properties of 
the solutions, we reconsider the above three cases, 
namely the cases where P(t) — 



4c' 



Pit) 



-y/f c tan(2 v / ^t), and fi(t) = tanh(2v^cSt) 
which correspond to the respective scattering lengths 
v(t) — v§ exp(— 2yJ— ca £), v(t) = vo/\cos(2y/cat)\, and 
v(t) = vo/ cosh(2v / — cat). The solution of equation (j4]) 
in each of these cases, for cs — and C3 7^ 0, respectively, 
may read 



i/>(z, t) = A s(t) sech[v^(C - Co)] e 



i(/3z 2 +fc£+ft) 



^OM) = A 5(t) tanh[-v^(C - Ci)]e : 



i(/3z 2 +/cz+fi) 



(21) 



(22) 



The time-dependent amplitude coefficient, s(t) = 
tit)~ x l 2 , depends on the typical forms of /3(t) into con- 
sideration. We have A = ^^fp 2 and Ci = Co - ] j=- 

Additionally, we recall that C = pz + q and Co = 
a — m(4c2|c4|). Without loss of generality, we set 

qi = 0. C2 
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A. The vanishing matter waves 




When j3 = ^ — ^, the solution of equation (jl]) is given 
in Eqs. ([2T]) and ([22]) . where p, /3, A: and Vt are time- 
dependent functions given in Eq. ([TO]) for the solution in 
Eq. ((2H), and in Eq. (dU) for the solution in Eq. ([22j) . As 
already said, these solutions represent bright and dark 
solitons, respectively. However they have common dy- 
namical behavior. The width of each of them is propor- 
tional to (pqC2)~ 1 / 2 e 4c/3 *, while the height is proportional 
to A = Aoe _2c ^. In this case, the matter waves have 
broadening and vanishing properties. As a matter of fact, 
with time, the height of the matter wave decreases while 
its width increases. However the numb er of atoms in the 

condensate, i.e. N = J \t/j\ 2 dz = — 4c ^° C2 -, remains un- 
changed during the propagation of the wave. A plot of 
this dynamical behavior is given in Figs. QJa) and [2^a) 
through the space-time evolution of the square magni- 
tude of the wave function. 

The kinematics of the wave can be obtained from 
t) = 0. Hence the motion of the center of mass, 
taken as the position that corresponds to the peak, is 
determined by the following equation: 



(23) 



for both dark and bright solitons. At longer times, the 
center of mass of the wave is driven towards the point 
z = — ^ = Zoo which corresponds to the effective trap 
center. In fact, when the gravitational field is considered, 
the minimum of the potential is no more on the magnetic 
trap axis z = 0, it moves to z^. Hence the gravitational 
field drives the wave from the center of the magnetic 
trap to a region around z^, where the wave should be 
confined. 

The velocity of the wave packet is z = ck\ e ~ 2 V-cat^ 
Hence k\ appears to be a measurement of the initial ve- 
locity of the wave. The velocity exponentially decreases 
with time. So the choice of the parameters A and k\ can 
seriously affect the dynamics of the matter waves, de- 
nouncing the role of the gravitational field in our analysis. 
The wave packet behaves like a static classical particle for 
k\ = and like a moving one for k\ ^ 0. A similar result 
was obtained in [25| within a special case where the grav- 
itational field is absent (A = 0) and C3 = 0. The accel- 
eration of the wave packet is z = —2^J—(?otk\ e ~ 2 V-cat^ 
The acceleration exponentially decreases with time. 

When /3 = sf—j^ tanh(2 v / — ca i), the solution of equa- 
tion (jlj) is given in Eqs. ([2T]) and ([22]) . where £, p, q, /?, 
k and £1 are time-dependent functions given in equation 
([T2]) for the solution in Eq. ([2T]) , and in equation ([T5]) for 
the solution in Eq. ([22]) . These solutions also represent 

growing matter waves. The height of each of these waves 

/ 2c 2 C \l/2 
is proportional to ( - ztztttt^^ 1 , and the width is 



vo cosh(2y / cat) J 

proportional to ^ cosh(2 v / ca t). The width of the soli- 
ton shortens exponentially with time while its height ex- 



ponentially grows. A display of this dynamical behavior 
can be found in Figs. [U^b) and I2fb) where we plot the 
space-time evolution of the wave in the system. The mo- 
tion of the center of mass of the matter wave is defined 
by the equation: 



z = —ki 



2a' 



(24) 



At longer times, the center of mass of the wave is driven 
towards the point z = ^ = z^ which corresponds to the 
effective trap center. The velocity of the wave packet is 
z = 2ck\ e - 2 V-cat^ wn j c h i s twice the velocity in the pre- 
vious case. We have observed that the solutions in the 
cases f3(t) = V^lf and f3(t) \/^t c tanh(2 v / ^5t), 
both corresponding to an expulsive trapping potential, 
present similar asymptotic behavior in time. In fact, 
when t — > 00 the width, height, and trajectory in both 
cases become identical. 



B. The growing matter waves 

When f3 = —\f^- c tan(2 v / cat), the solution of equa- 
tion ((4]) is given in Eqs. ([2T]) and ([22]) . where i, p, q, /?, 
k and ft are time-dependent functions given in equation 
([TT]) for the solution in Eq. ([2T]) . and in equation ([T7]) 
for the solution in Eq. ([22]) . The corresponding solutions 
represent growing matter waves. The height of each of 

these matter waves is proportional to ( VQ | C os(2°v^a t) \ ) ' 
and the width is proportional to cos(2^/cat)\. In the 
safe time interval, the matter wave becomes thinner and 
higher. We portray in Figs. [TJc) and[2fc) this dynamical 
behavior. Close to the cut-off time which is t\ = 15.7, the 
exponential increase in the amplitude of the wave is so 
strong that a "collapse" of the wave may occur. However, 
by changing the expression of the parameter j3(t) (which 
amounts to changing the expression of the s-wave scat- 
tering length) in a small time interval that contains each 
singular time t n = ^ 2 ^^ 7r , the propagation of the mat- 
ter wave can be kept. It can be changed to j3(t) = ^ — 1| 
or /3(t) = y/— H tanh(2 v 7 — cat). In this case, the width 
and the height of the wave oscillate in time. Figures [3ja) 
and [3](b) show the long-time propagation of the matter 
waves in this case. 

The motion of the center of mass of the matter wave 
is defined by the equation: 



c . / 1 — \ A 

sm(2v cat) . 

a 2a 



(25) 



The velocity of the wave packet is z = 2ck\ cos(2^/cat) 
while its acceleration is z = 4Vc 3 ak± cos{2yJcat). This 
means that the wave oscillates in time with frequency 
/ = also equivalent to (£ n +i — ^n-i) -1 - These os- 
cillations are performed around the position z = — 
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X (b) 

FIG. 3: (Color online) Plot of the long-time evolution of the 
matter waves in the case where /3(t) = — tan(2y / ca t) 
for (a) the solution in Eq. ([T3|h and (b) the solution in Eq. 
([T9]) . We use f3(t) — y / — ^ tanh(2^/— cat) in small time 

intervals that contain each singular time t n = , with n 

being any integer. The parameters are a = ±0.005, po = 1.0, 
— —1.0, a = 0, c = 0.5, C2 = 1, C4 = i/o, &i = 0.8, gi = 0, 
A = 0.01, and 6 = 1. All quantities are dimensionless. 

which corresponds to the effective trap center set by the 
gravitational field. 

In comparison with the results obtained in [4] , we see 
that the role of gravity would not be the same, even qual- 
itatively, for different traps of the BEC system. In fact, 
the bias magnetic field may be the analog of gravitational 
field since both fields are represented by the linear term 
in the trapping potential. From the results of [4|, we 
infer that the kink solitons created in a bias magnetic 
field alone behave like a classical particle in a pure free 
fall motion led by the "gravity" in the (z,t) space. In 
the present case, the gravitational field is not alone. The 
presence of the parabolic magnetic potential changes the 
effect of the gravitational field. For instance, when the 
initial speed of the bright or dark soliton is zero, then 
both the velocity and the acceleration at any time are 
zero too. 

The present study suggests three ways to generate 
bright and dark solitons in BEC systems by time- 
varying the s-wave scattering length (through the Fes- 
hbach resonance) without changing the magnetic poten- 



tial. This can be done by tuning the scattering length 
to go/\cos(2y/cat)\ when the condensate is confined in 
an attractive magnetic trap, i.e. a is positive. We may 
also tune the scattering length to go/ cosh(2 v /— ca i), or 
go ex.p(—2^—ca i), in the case of expulsive magnetic trap, 
i.e. a is negative. 

V. CONCLUSION 

In conclusion, we have considered the GP equation 
with time-dependent cubic nonlinearity which describes 
the dynamics of the BEC matter-waves in a magnetic 
field and under the effect of a homogeneous gravita- 
tional field. With the help of the extended tanh-function 
method, we have obtained a solution which has as special 
cases the bright and dark solitons. As has been discussed, 
these solitons can be generated by properly tuning the 
5-wave scattering length of the condensed particles, de- 
pending on whether the magnetic trapping is attractive 
or expulsive. The dynamics and kinematics of these mat- 
ter waves have been presented and discussed. 

We have found that the gravity reshapes the repel force 
of the magnetic trap, and then drives the matter waves 
towards the region around the position z = — The 
matter waves may remain in that position for scatter- 
ing lengths go/ cosh(2v"— ca t) or go exp(— 2\J— ca t), and 
may oscillate around it for go/\ cos(2y / ca t)\. By compar- 
ing the results obtained here with those of [4], we have 
found that the role of gravitational field depends on the 
type of trap in which the condensate is confined. 

The study of the dynamics and stability of a BEC un- 
der the effect of very strong gravitational field, that could 
occur (in a speculative way) for instance close to black 
holes, appears to pose an interesting issue to investigate 
in future works. 



Acknowledgments 

Part of this work has been done during the Short 
Term Visit of EW within the CMSP Section of the Ab- 
dus Salam ICTP (Italy). EW acknowledges the support 
from the Government of India, through the CV Raman 
International Fellowship for African Researchers. HV 
acknowledges the support from CNPq (Brazil) and the 
DFG Research Training group 1620 "Models of Gravity" . 
AM Thanks the Abdus Salam ICTP for financial support 
through the Associateship program. KP acknowledges 
DST, DAE-BRNS, UGC, the Government of India, for 
financial support through major projects. 



[1] V. E. Zakharov and S. V. Nazarenko, Physica D 201, U. K., 1995). 

203 (2005). [3] A. Mohamadou, E. Wamba, D. Lissouck, and T. C. Ko- 
[2] A. Griffin, D. W. Snoke, and S. Stringari, Bose-Einstein fane, Phys. Rev. E 85, 046605 (2012). 

Condensation (Cambridge University Press, Cambridge, [4] E. Wamba, T. C. Kofane, and A. Mohamadou, Chin. 



8 



Phys. B 21, 070504 (2012). 
[5] L.-C. Zhao, Z.-Y. Yang, T. Zhang, and K.-J. Shi, Chin. 

Phys. Lett. 26, 120301 (2009). 
[6] Z. Yan, K. W. Chow, and B. A. Malomed, Chaos, Soli- 
tons and Fractals 42, 3013 (2009); U. Al Khawaja, J. 

Math. Phys. 51, 053506 (2010). 
[7] R. Murali and K. Porsezian, Physica D 239, 1 (2010). 
[8] Y. Kawaguchi, M. Nakahara, and T. Ohmi, Phys. Rev. 

A 70, 043605 (2004). 
[9] A. Mohamadou, E. Wamba, S. Y. Doka, T. B. Ekogo, 

and T. C. Kofane, Phys. Rev. A 84, 023602 (2011). 
[10] J. I. Rivas and A. Camacho, Mod. Phys. Lett. A 26, 481 

(2011). 

[11] G. Bertone, D. Hooper, and J. Silk, Physics Reports 405, 
279 (2005). 

[12] H. Velten and E. Wamba, Phys. Lett. B, 709 1 (2012); 

S. J. Si, Phys. Rev. D 50, 3650 (1994); C. G. Bohmer 

and T. Harko, JCAP 0706, 025 (2007). 
[13] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86, 2353 

(2001). 

[14] F. Dalfovo, S. Giorgini, L.P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[15] F. Kh. Abdullaev, B. B. Baizakov, S. A. Darmanyan, V. 

V. Konotop, and M. Salerno, Phys. Rev. A 64, 043606 

(2001). 



[16] P. S. Julienne, F. H. Mies, E. Tiesinga, and C. J. 

Williams, Phys. Rev. Lett. 70, 1880 (1997). 
[17] A. Camacho, L. F. Barragan-Gil, and A. Macfas, Cent. 

Eur. J. Phys. 8, 717 (2010). 
[18] A. J. Legget, Rev. Mod. Phys. 73, 307 (2001). 
[19] J. R. Ensher, PhD thesis, University of Colorado (Boul- 
der, USA, 1998). 
[20] I. K. Kulikov, |arXiv:cond-mat/0205330^ 3 [cond- 

mat.stat-mech] (2002); I. K. Kulikov, Int. J. Theor. 

Phys. 41, 1481 (2002). 
[21] B. P. Anderson and M. A. Kasevich, Science 282, 1686 

(1998); W. Zhang and D. F. Walls, Phys. Rev. A 57, 

1248 (1998). 

[22] A.-X. Zhang and J.-K. Kue, Phys. Rev. A 75, 013624 
(2007); F. Kh. Abdullaev, R. M. Galimzyanov, and Kh. 
N. Ismatullaev, J. Phys. B: At. Mol. Opt. Phys. 41, 
015301 (2008). 

[23] Ph. Courteille, R. S. Freeland, D. J. Heinzen, F. A. van 
Abeelen, and B. J. Verhaar, Phys. Rev. Lett. 81, 69 
(1998); E. Timmermans, P. Tommasini, M. Hussein, and 
A. Kerman, Physics Reports 315, 199 (1999). 

[24] G. Nandi, R. Walser, E. Kajari, and W. P. Schleich, Phys. 
Rev. A 76, 063617 (2007). 

[25] H.-M. Li, Chin. Phys. 15, 2216 (2006). 



